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Abstract. This article concerns the computational problem of count- 
ing the lattice points inside convex polytopes, when each point must 
be counted with a weight associated to it. We describe an efficient al- 
gorithm for computing the highest degree coefficients of the weighted 
Ehrhart quasi-polynomial for a rational simple polytope in varying di- 
mension, when the weights of the lattice points are given by a polynomial 
function h. Our technique is based on a refinement of an algorithm of 
A. Barvinok in the unweighted case (i.e., h = 1). In contrast to Barvi- 
nok's method, our method is local, obtains an approximation on the 
level of generating functions, handles the general weighted case, and 
provides the coefficients in closed form as step polynomials of the di- 
lation. To demonstrate the practicality of our approach we report on 
computational experiments which show even our simple implementation 
can compete with state of the art software. 
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1. Introduction 

Computations with lattice points in convex polyhedra arise in various 
areas of computer science, mathematics, and statistics (see e.g., [I3l I23] [32] 
and the many references therein). Given p, a rational convex polytope in M. d , 
and h(x), a polynomial function on R" (often called a weight function) , this 
article considers the important computational problem of computing, or 
estimating, the sum of the values of h(x) over the lattice points belonging 
to p, namely 

s( P ,h)= 

The function S(p,h) has already been studied extensively in the unweighted 
case, i.e., when h{x) takes only the constant value 1 (in that case of course 
S(jp, 1) is just the number of lattice points of p). Many papers and books 
have been written about the structure of that function (see, e.g., [TTJ [13] 
and the many references therein). Nevertheless, in many applications h(x) 
can be a much more complicated function. Important examples of such a 
situation appear, for instance, in enumerative combinatorics PQ, statistics 
|27[I22|. symbolic integration [2] and non-linear optimization |24] , Still, only 
a small number of algorithmic results exist about the case of an arbitrary 
polynomial h. 

It is well-known that when the polyhedron p is dilated by an integer 
factor n € N, we obtain a function of n, the so-called weighted Ehrhart 
quasi-polynomial of the pair (p, h), namely 

d+M 

S(np,h) = ^] h(x) = ^] E m (n mod q) n m . 

x£npnZ d rn=0 

This is a quasi-polynomial in the sense that the function is a sum of mono- 
mials up to degree d + M, where M = deg/i, but whose coefficients E m are 
periodic functions of n. The coefficient functions E m are periodic functions 
with period q, where q E N is the smallest positive integer such that qp 
is a lattice polytope, i.e., its vertices are lattice points. We will make this 
more precise later (we recommend [11] [13] for excellent introductions to this 
topic). 

To begin realizing the richness of S(np,h), note that its leading highest 
degree coefficient E^+m (which actually does not depend on n) is precisely 
equal to f h(x)dx, i.e., the integral of h over the polytope p, when h is 
homogeneous of degree M. These integrals were studied in [6], [7] and more 
recently at pj. Still most other coefficients are difficult to understand, even 
for easy polytopes, such as simplices (see [21] for a survey of results and 
challenges). The key aim of this article is to achieve the fast computation 
of the first few top-degree (weighted) coefficients E m via an approximation 
of S(np, h) by a quasi-polynomial that shares the highest coefficients with 
S(np, h). 



HIGHEST EHRHART COEFFICIENTS 



3 



We now explain the known results achieved so far in the literature. It 
is important to stress that computing all the coefficients E m for m = 
0, . . . , d + M is an NP-hard problem, thus the best one can hope for theoret- 
ical results is to obtain an approximation, as we propose to do here. Until 
now most results dealt only with the unweighted case, i.e., h{x) = 1 and 
we summarize them here: A. Barvinok first obtained for lattice polytopes p 
a polynomial-time algorithm that for a fixed integer fco can compute the 
highest ko coefficients E m (see [8]). For this he used Morelli's identities [33] 
and relied on an oracle that computes the volumes of faces. 

Later, in [TO], Barvinok obtained a formula relating the k highest degree 
coefficients of the (unweighted) Ehrhart quasi-polynomial of a rational poly- 
tope to volumes of sections of the polytope by certain affine lattice subspaces 
of codimension < k. As a consequence, he proved that the k highest degree 
coefficients of the unweighted Ehrhart quasi-polynomial of a rational simplex 
can be computed by a polynomial algorithm, when the dimension d is part 
of the input, but k is fixed. More precisely, given a dilation class n mod q 
with n € N, Barvinok's algorithm computes the numbers E m (n mod q) by 
an interpolation technique. However, neither a closed formula for these 
E m , depending on n, nor a generating function for the coefficients became 
available from [TO]. In fact, in that article [10\ Section 8.2] the question of 
efficiently computing such a closed form expression was raised. 

A key point of both Barvinok's and our method is the following. The sum 
S(p,h) has natural generalizations, the intermediate sums S L (p,h), where 
L C V = M. d is a rational vector subspace. For a polytope p C V and a 
polynomial h(x) 

S L (p,h) = Y, [ Ky)dy, 
x Jpn(x+L) 

where the summation index x runs over the projected lattice in V/L. In 
other words, the polytope p is sliced along lattice affine subspaces parallel 
to L and the integrals of h over the slices are added up. For L = V, there 
is only one term and S v (p, h) is just the integral of h(x) over p, while, for 
L = {0}, we recover S(p,h). Barvinok's method in [TO] was to introduce 
particular linear combinations of the intermediate sums, 

^KL)S L (p,h). 

It is natural to replace the polynomial weight h{x) with an exponential 
function x i— > e^' x \ and consider the corresponding holomorphic functions 
of ^ in the dual V*. Moreover, one can allow p to be unbounded, then the 
sums 

s L (pm=J2 I e{u)d y 

x Jpn(x+L) 

still make sense as meromorphic functions on V*. The map p i— > S L (p)(S,) is 
a valuation. 
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In [14], it was proved that a version of Barvinok's construction on the level 
of generating functions, namely X^Lg£ ^(L)S L (p)(£), approximates 5(p)(^) 
in a certain ring of meromorphic functions (a precise statement is given 
below). The proof in [T4] relied on the Euler-Maclaurin expansion of these 
functions. Another proof, using the Poisson summation formula, will appear 
in p2]. 

In the present article, we introduce a simplified way to approximate 
S(p)(t;) for the case of a simplicial afhne cone p = s + c, which levels the 
way for a practical and efficient implementation. Via Brion's theorem, it 
is sufficient to sum up these local contributions of the tangent cones of the 
vertices. We present a method for computing the highest degree coefficients 
of the Ehrhart quasi-polynomial of a rational simple polytope, by applying 
the approximation theorem to each of the cones at vertices of p. The com- 
plexity depends on the number of vertices of the polytope, and thus if the 
simple polytope is presented by its vertices (rather than by linear inequali- 
ties), we obtain a polynomial-time algorithm. In particular, the algorithm 
is polynomial-time for the case of a simplex. We obtain the Ehrhart co- 
efficient functions E m (n mod q) in a closed form as step polynomials, i.e., 
polynomials in n whose coefficients are modular expressions (Ci n ) mod qi 
with integers d and qi, for example (2n) mod 3. Having a closed formula 
available considerably strengthens Barvinok's result in [10] even in the un- 
weighted case h = 1. 

The structure of this paper is as follows. In Section [2] we first present 
some necessary preliminaries. Section [3] explains the intermediate generat- 
ing function S L (p)(£) in more detail. Then we show how to use a grading 
of S(p)(!;) to extract the highest degree coefficients of the weighted Ehrhart 
polynomial in the case of a lattice polytope. This motivates the approxima- 
tion results for generating functions. In Section [4] we give a simple proof of 
the approximation theorem of [14] , in the case of a simplicial cone (see The- 
orem l24l) . The theorem uses the notion of a patching function (essentially 
a form of Mobius inversion formulas described in Subsection 14. ip . We ex- 
hibit an explicit and easily computable such patching function. Using these 
tools, we show in Section [5] that the approximation for a cone s + c (on the 
level of generating functions) can be computed efficiently as a closed for- 
mula. The formula makes the periodic dependence on the vertex s explicit. 
Finally, in Section [6] we give the polynomial-time algorithm to compute 
the coefficients E m (n mod q) as step polynomials. Our main result (Theo- 
rem [37]) says that, for every fixed number ko, there exists a polynomial-time 
algorithm that, given a simple polytope p of arbitrary dimension, a lin- 
ear form £ G V* , a nonnegative integer M, computes the highest ko + 1 
coefficients Em +d-ko >•••> ^M+d of the weighted Ehrhart quasi-polynomial 
S(np, h = i M ) in the form of step polynomials. 

Four comments are in order about the applicability and potential practi- 
cality of the main results: First, although the weight h used in Theorem [37] 
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is a power of a linear form, as is carefully explained in [2], one can obtain 
similar complexity of computation for polynomials that depend on a fixed 
number of variables, or with fixed degree (Corollary 04]). Second, it is also 
worth noting that using perturbations (see e.g., [21]), triangulations [26] or 
simplicial cone decompositions of polyhedra (see, e.g., [30J), one can extend 
computations from simple polytopes to arbitrary polytopes. Third, since 
our approximation is done at the level of generating functions, it extends 
the complexity result from |10j to the weighted case. Finally, at the end of 
the article we report on experiments using a simple implementation of the 
algorithm in Maple, demonstrating it is competitive with more sophisticated 
software tools. This indicates a potential to use this algorithm for experi- 
mentally verifying conjectures on the positivity of the Ehrhart coefficients of 
certain polytopes, for examples where the computation of the full Ehrhart 
polynomials is out of reach. The algorithms presented here require a rich 
mixture of computational geometry and algebraic-symbolic computation. 



2. Preliminaries 

2.1. Rational convex polyhedra. We consider a rational vector space V 
of dimension d, that is to say a finite dimensional real vector space with 
a lattice denoted by A. We will need to consider subspaces and quotient 
spaces of V, this is why we cannot simply let V = K rf and A = A point 
v £ V is called rational if there exists a non-zero integer q such that qv G A. 
The set of rational points in V is denoted by Vq. A subspace L of V is called 
rational if L n A is a lattice in L. If L is a rational subspace, the image of 
A in y/i is a lattice in V/L, so that V/L is a rational vector space. The 
image of A in V/L is called the projected lattice. A rational space V, with 
lattice A, has a canonical Lebesgue measure dx = dm\(x), for which V/A 
has measure 1. 

A convex rational polyhedron p in V (we will simply say polyhedron) is, by 
definition, the intersection of a finite number of closed half spaces bounded 
by rational affine hyperplanes. We say that p is full- dimensional (in V) if 
the affine span of p is V. 

In this article, a cone is a polyhedral cone (with vertex 0) and an affine 
cone is a translated set s + c of a cone c. A cone c is called simplicial if 
it is generated by independent elements of V. A simplicial cone c is called 
unimodular if it is generated by independent integral vectors v\, . . . ,Vk such 
that {m, . . . , Vk} can be completed to an integral basis of V. An affine cone a 
is called simplicial (respectively, simplicial unimodular) if the associated 
cone is. A polytope p is a compact polyhedron. The set of vertices of p is 
denoted by V(p). For each vertex s, the cone of feasible directions at s is 
denoted by c s . For details in all these notions see, e.g., [TT] . 
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2.2. Generating functions: Exponential sums and integrals. 

Definition 1. We denote by H(V*) the ring of holomorphic functions de- 
fined around G V* . We denote by A4(V*) the ring of meromorphic func- 
tions defined around G V* and by A4e(V*) C A4(V*) the subring consist- 
ing of meromorphic functions 0(£) which can be written as a quotient of a 
holomorphic function and a product of linear forms. 

This paper relies on the study of important examples of functions in 
M.e(V*), the following continuous and discrete generating functions I(p, A) 
and 5(p, A) associated to a convex polyhedron p. Both have an important 
additivity property which makes them valuations (see Chapter 8] or the 
survey [12] for a detailed presentation, here we summarize the essentials). 

Definition 2. Let M be a vector space. A valuation F is a map from the set 
of polyhedra p C V to the vector space M such that whenever the indicator 
functions [pi] of a family of polyhedra pi satisfy a linear relation r i [pi] 
. then the elements F(pi) satisfy the same relation Yli r i^{Pi) = 0- 

Proposition 3. There exists a unique valuation /(■, A) which associates to 
every polyhedron p C V a meromorphic function I(p, A) G ftA(.(V*), so that 
the following properties hold: 

(i) // the polyhedron p is not full- dimensional or if p contains a straight 
line, then I(p, A) = 0. 

(ii) If^&V* is such that e^' x ^ is integrable overp, then 

I(p,A)(0= fe^dm A (x). 
Jp 

(iii) For every point s G Vq, one has 

/( S + p,A)(0 = e^>/(p,A)(0. 

We will call I(p,A)(£) the continuous generating function of p. 

Proposition 4. There exists a unique valuation ^(^A) which associates to 
every polyhedron p C V a meromorphic function S(p, A) G M.i(V*), so that 
the following properties hold: 

(i) If p contains a straight line, then S(p, A) = 0. 

(ii) If ^ G V* is such that e^' x ^ is summable over the set of lattice points 
ofp, then 

5(p,A)(0= £ e<^>. 
xe pnA 

(iii) For every point s G A, one has 

S(s + p,A)(0 = e^S(p,A)(0- 
5"(p,A)(^) is called the (discrete) generating function ofp. 
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2.3. Brion's theorem. A consequence of the valuation property is the fol- 
lowing fundamental theorem. It follows from the Brion-Lawrence-Varchenko 
decomposition of a polyhedron into the supporting cones at its vertices 
[18^ 111]: see also [19] . Proposition 3.1, for a more general Brianchon-Gram 
type identity. 

Theorem 5. Let p be a polyhedron with set of vertices V(p). For each 
vertex s, let c s be the cone of feasible directions at s. Then 

S(p,A)= S(s + c s ,A). 

sev{p) 

2.4. Notations and basic facts in the case of a simplicial cone. For all 

of the notions below see [TT]. Let Uj £ A, i = 1, . . . , d be linearly independent 
integral vectors and let c = X^j=i ^+ v i be the cone they span. 

Definition 6. The fundamental parallelepiped b of the cone (with respect 
to the generators V{, i = 1, . . . , d) is the set 

d 

1=1 

Note that the set has a half-open boundary. We immediately have: 
Lemma 7. Let s £ V . Then 

/(.+., axo-.m ^^w , (1) 
n<=i (^ v i) 

where vo1a(£>) is the volume of the fundamental parallelepiped with respect 
to the Lebesgue measure dm\ defined by the lattice. 

If V = M d and A = Z d , then vol A (b) = |det(ui, . . .,v d )\, and so 

nii<Oi> 

We also recall the following elementary but crucial lemma. 

Lemma 8. (i) The affine cone (s + c) n A is the disjoint union of the 

translated parallelepipeds s + b + v, for v G Ylj=i ^Vj. 
(ii) The set of lattice points in the affine cone s + c is the disjoint union of 

the sets x + Yli=i when x runs over the set (s + b) n A. 
(hi) The number of lattice points in the parallelepiped s + b is equal to the 

volume of the parallelepiped with respect to the Lebesgue measure dm^ 

defined by the lattice, that is 

Card((s + b) n A) = vol A (b). 
In particular, when V = M. d and A = then 

Card((s + b) HZ d ) = |det(«i,...,u rf )|. 
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The study of the generating function S(s + c, A)(£) of the affine cone s + c 
will be a crucial tool. It relies on expressing S(s + c, A)(£) in terms of the 
generating function S(s + b, A)(£) = Ylxe(s+b) nA &^' x ' of the fundamental 
parallelepiped. Lemma [8] (ii) immediately gives: 

Lemma 9. 

S(s + c, A)(0 = S(s + b, A)(£) , 1 (3) 

ni=i(i-e^^) 

Example 10. Consider the case where V = M and A = Z. Let c = K + . Let s G K, 
then the "fractional part" {s} G [0, 1[ is defined as the unique real number such 
that s — {s} G Z. Then the unique integer s in s + b is s + {— s}, and so equations 
O and © give 

J(a + c,Z)(0 =e« s ^: an d 5( s + c , Z)(0 = c«( s+ ^ s » — L^. 

^ re 

3. Key ideas of the approximation theory 

3.1. Weighted Ehrhart quasi-polynomials. Let p C V be a rational 
polytope and let be a polynomial function of degree M on V. We 

consider the following weighted sum over the set of lattice points of p, 

zepnA 

When p is dilated by a non-negative integer n G N, we obtain the weighted 
Ehrhart quasi-polynomial of the pair (p,h). 

Definition 11. Xei g be the smallest positive integer such that qp is a lattice 
polytope. The we define the Ehrhart quasi-polynomial E(p,h;n) and its 
coefficients E m (p, h; n mod q) by 

d+M 

E(p,h;n)= H x ) = E E m(p, h; n mod q) n m . 

xGrapnA m=0 

We note that the coefficients E m depend on n, but they actually depend 
only on n mod q, where q is the smallest positive integer such that qp is a 
lattice polytope. If h(x) is homogeneous of degree M, the highest degree 
coefficient E^+m is equal to the integral J h{x) dx (see [2] and references 
therein) . 

We concentrate on the special case where the polynomial h(x) is a power 
of a linear form 

h(x) = ^ 
h[X) Ml • 

This is not a restriction because any polynomial can be written as a linear 
combination of powers of linear forms. In fact, as discussed in [2j, whenever 
the polynomial h(x) is either of fixed degree or only depends on a fixed 
number of variables (possibly after a linear change of variables), then only a 
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Figure 1. The example triangle t and its dilation 6t 



polynomial number of powers of linear forms are needed, and such a decom- 
position can be computed in polynomial time. We introduce the following 
notation. 

Definition 12. Let q be as above. We define the Ehrhart quasi-polynomial 
E(p, £, M; n) and the coefficients E m (p, £, M; n mod q) for m = 0, . . . , M + d 
by 



E(p,£,M;n)= £ 



Ml 



E m (p,£,M;nmod 



i n 



xCnpnA m=0 

It will be convenient in this paper to introduce the following notations. 
For a positive integer q € N and a real number n 6 R, we write 

[njg := g[|nj 6 gZ, {n} 9 := (ra mod g) G [0,g), 

which give the unique decomposition 

n = N<? + i n }q- 

By \n\ := [nji and {n} := we obtain the ordinary "floor" and "frac- 

tional part" notations. Finally, |~n] := — [_— n\ is the "ceiling" notation. 

Example 13. Consider the rational triangle t with vertices (0,0), (^,0), and 
(m> Ti) as shown in Figure [TJ Let us compute E(i, £, M; n) for this small example. 
Note that the integer q such that qt is a lattice polytope is q = 28. 

In what follows consider powers of the linear form £ = x + y as weights for the 
lattice points. When the power M = 0, then we obtain a constant weight and the 
quasi-polynomial E(i, £, 0; n) counts the lattice points inside the various dilations 
of t: it is given by the following formula: 



25 2 ( 5 , , 5 
784 n + l-392 {5n}28+ 14 1 " 



l + ^({5n} 28 f-^{5n} 28 



Indeed, when n = 1 (no dilation) there is only one lattice point and the formula 
above reduces to 



1089 33 1 2 

"78T " 392 |5}28 + 784 ({5|28) 



-__--_-5+J-25=l 
784 392 784 
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When we dilate the same triangle six times, i.e., n = 6, we obtain four lattice 
points. 

841 2Q 1 841 2Q 1 

Next let us take M = 1, in that case the lattice point (a, b) is counted with weight 
a + b. In this case the top coefficient is equal to the integral of the linear form 
£ = x + y over t. The quasi-polynomial E(t, £, l;n) is given by the following 
formula: 

125 3 / 25 75 \ 2 J 5 /rr _ 1 , 2 15 rc _ 1 , 25 



K,4«4 n + l-5488 {5n}28 + 784j " + (,5488 ({5 nU) ~ 392 {5 nU + 84 

1 3 5 

({5n}28) +3 5 7({5n} 2 s) -^7{5n} 2 8 



16464 Vl ' ' 784 VL J ' 84 

Substitute again n — 1 in the expression, to obtain 

275 1685 13 2 1 3 n 

686 " 16464 {5}a8 + 2744 ({5}2s) " 16464 ({5}2s) = °' 

Note that since only the lattice point (0,0) lies within the triangle at n = 1, the 
quasipolynomial must evaluate to zero. 

In practice, it is impossible to compute E(p,£,M;n) except when p is of 
small dimension (and M relatively small). Thus we restrict our ambitions: 

Let us fix a number ko. Our goal will be to compute the ko + 1 highest 
degree coefficients E m (p, £, M; n mod q), for m = M + d, . . . , M + d — ko- 
We will be able to give a polynomial time algorithm to do so. 



3.2. Grading of the generating functions. A key property that we will 
make use of in this article is the following grading of Mi: A function </>(£) G 
AiiiV*) has a unique expansion into homogeneous rational functions 

m>mo 

where the summands <p[ m ]{i) have degree m as we define now: If P is a 
homogeneous polynomial on V* of degree p, and D a product of r linear 
forms, then € is an element in Aie(V*) homogeneous of degree m = p — r. 
For instance, g is homogeneous of degree 0. On this example we observe 
that a function in Mi(V*) which has non- negative degree terms need not 
be analytic. 

In particular, consider the generating function S(s + c, A)(£) for a simpli- 
cial cone c. By Lemma [U 

S(s + c, A)(0 = S(s + b, W ^^^y 

Thus, S(s + c, A) € Aii(V*), and so it admits a decomposition into homo- 
geneous components: 
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Lemma 14. 

S(s + c, A) (0 = 5(s + c, A) M] (0 + S(s + c, A) (0 + • • • , (4) 

and the lowest degree term S(s + c, A)j_ d ](£) is equal to J(c, A)(£), i.e., the 
integral over the unshifted cone c. 

Proof. We write 

rr 1 -rr & Vj) 1 (5) 



The function yr^F is holomorphic with value — 1 for s = 0. Thus we have 
S(s + c, A) € Me(V*). The value at £ = of the sum over the parallelepiped 
is the number of lattice points of the parallelepiped, that is voL\(b)- This 
proves the last assertion. □ 

3.3. Sketch of the method for lattice polytopes. We will now explain 
the key point of our method, with the simplifying assumption that the ver- 
tices of the polytope are lattice points. We will show that the highest degree 
coefficients of the weighted Ehrhart polynomial can be read out from an ap- 
proximation of the generating functions of the cones at vertices. In Section [J] 
we will study this approximation, and in Section [5] we will show how to effi- 
ciently compute it. Then, in Section [6j we will come back to the computation 
of Ehrhart coefficients for the general case of rational polytopes. 



Proposition 15. Let p be a lattice polytope. Then, for k > 0, we have 



{£,s) M+d - k 



E M+d - k (p,Z,M)= ^ M + d-kV S ^- d+ ^°- (6) 

sGV(p) V 



The highest degree coefficient is just the integral 



f (£ x) M 
E M+d (p,Z,M) = J Ai ^-dx. 



Remark 16. As functions of£, the coefficients E m (p, £, M) are polynomial, 
homogeneous of degree M. However, in ([6]), they are expressed as linear 
combinations of rational functions of £, whose poles cancel out. 

Proof of Proposition \lb~[ The starting point is Brion's formula. As the ver- 
tices are lattice points, we have 

E ^ ,X) = E S(s + c s )(0= E e«-'>S(c.)(0. (7) 
xepnA seV(p) sev(p) 

When p is replaced with np, the vertex s is replaced with ns but the cone 
c s does not change. We obtain 

E e< ^ } = E ^ s) s(c s )(o- 

xenpnA sGV(p) 
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We replace £ with i£, 

£ e 4 ^ = £ e ntiM S(CsM)- 
xenpnA seV(p) 

The decomposition into homogeneous components gives 

S(tsM) = t- d I(c s m + t- d+1 S(c s ) hd+1] (0 + ■■■ + t k S(c s ) [k] (0 + • • • . 
Hence, the i^-term in the right-hand side of the above equation is equal to 

M+d It „\M+d-k 

Thus we have 



V - V n M+d ^ ,S ^ M+d I(c )(f) 

xenpnA seV(p) 



+ n M+d - 



1 (Z,s) M+d - 1 



(M + d-l)! S(Cs) [- rf+1 ' (0 + ' ' ' + *kWfl- ( § ) 

^ /t X \M 

From this relation, we read immediately that }2xenp nA aii ^ s a polynomial 
function of n of degree M + d, and that the coefficient of n M+d ~ k is given 
by ©. The highest degree coefficient is given by 

It „\M+d 

sG v( P ) v ^ 

Applying Brion's formula for the integral, this is equal to the term of £- 
degree M in I(p)(£), which is indeed the integral f ^'m\ ^ 

From Proposition [T5J we draw an important consequence: in order to 
compute the ko + 1 highest degree terms of the weighted Ehrhart polyno- 

/£ X \M 

mial for the weight h(x) = , we only need the ko + 1 lowest degree 

homogeneous terms of the meromorphic function S(c s )(£), for every vertex 
s of p. We compute such an approximation in Section [JJ it turns out to be 
sufficient also in the general case of a rational polytope. 

3.4. Intermediate generating functions. To obtain the approximation, 
we study generating functions which interpolate between the integral I(p, A) 
and the discrete sum S(p,A). This trend of ideas was first discussed by 
Barvinok in [10]. Let L be a rational subspace of V. To any polyhedron p we 
associate a meromorphic function S L (p,A)(£) G A4(V*), which is, roughly 
speaking, obtained by slicing p along afhne subpaces parallel to L through 
lattice points, and adding the integrals of e^'^ along the slices. Recall that 
the quotient space V/L is endowed with the projected lattice Ay/L- 
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Proposition 17. Let L (IV be a rational subspace. There exists a unique 
valuation S L (-,A) which to every rational polyhedron p C V associates a 
meromorphic function with rational coefficients S (p, A) € M(V*) so that 
the following properties hold: 

(i) If p contains a line, then 5 (p, A) = 0. 

(ii) 

S L (p,A)(0= £ / e^dy, (9) 

for every £ € V* such that the above sum converges. 
(iii) For every point s € A, we have 

S L (s + p,A)(0=e^S L (p,A)(0. 

We call the function S L (p,A) an intermediate generating function. The 
proof is entirely analogous to the case L = {0}, see Theorem 3.1 in [12] . and 
we omit it. 

For L = {0}, we recover the valuation S. For L = V, we have S (p, A) = 
I(p,A). In particular, if p is not full-dimensional, then S v (p, A) = 0. 

If p is compact, the meromorphic function S L (p, A)(£) is actually regular 
at £ = 0, and its value for £ = is the Q- valued valuation E L ± (p) considered 
by Barvinok [10] . 

Remark 18. The function S L (p, A) is actually an element of M.i(V*), just 
like the functions S{p, A) and I(p, A). This follows from an interesting de- 
composition that allows to write S L as a combination of terms using S and I 
for certain cones. This and other properties of the valuation S L (-,A) will be 
discussed in a forthcoming article [1] . 

4. Approximation of the generating function of a simplicial 

affine cone 

Let c C V be a simplicial cone with integral generators Vj, j = l,...,d, 
and let s E Vq. Let ko < d. In this section we will obtain an expression for 
the ko + 1 lowest degree homogeneous terms of the meromorphic function 
S(s + c)(£). Recall that if c is unimodular, the function S(s + c)(£) has a 
"short" expression 

S(» + c)(«=e«-->n T -I ST , 
i=i 

where Vi, i = 1, . . . , d are the primitive integral generators of the edges and 
s is the unique lattice point in the corresponding parallelepiped s + b. This 
is a particular case of Lemma [8] 

When c is not unimodular, it is not possible to compute efficiently the 
first ko terms of the Laurent expansion of the function S(s + c)(£), if ko 
is part of the input as well as the dimension d. In contrast, if ko is fixed, 
we are going to obtain an expression for the terms of degree < — d + ko 



14 V. BALDONI, N. BERLINE, J. A. DE LOERA, M. KOPPE, AND M. VERGNE 



which only involves a discrete summation over cones in dimension < ko and 
determinants. For example, the lowest degree term is |det(i>j)| f\j • 

4.1. Patching functions. For constructing the approximation, we will use 
a patching function. For / C {1, . . . , d}, we denote by Lj the linear span of 
the vectors Vi, i € I and by Lj- C V* the orthogonal subspace. We denote 
by I c the complement of I in {1, . . . , d}. 

Definition 19. We denote by J> do the set of subsets I C {l,...,d} of 

cardinality \I\ > do. A function I i-> X(I) on J> do is called a patching 
function if it satisfies the following condition. 

U ^1= E Ki)[Lf]. (10) 



IP 7 d " , - ■ 



Remark 20. The family of sub spaces Lj, \I\ > do is closed under sum, and 
the family of orthogonals Lj- is closed under intersection. The value X(I) 
plays the same as role as the Mobius function n{L) for Lj- that Barvinok 
[TO"! section 7] computes algorithmically for a certain family of subspaces L 
by walking the poset. From this discussion, it follows that patching func- 
tions do exist. The precise relation between Barvinok 's construction and the 
construction of the present paper will be studied in the forthcoming paper [5 J . 

We will compute a canonical patching function below, in Proposition [27J 
Let us state some interesting properties. 

Lemma 21. Let 1 1-> X(I) be a function on J> do - The following conditions 
are equivalent. 

(i) X is a patching function. 

(ii) E/ej-d JQ i A(-f) = 1 for every I € J~ 



>d - 



(hi) For 1 < i < d, let Fi{z) G be a formal power series (in one 

variable) with constant term equal to 1. Then 



mod terms of z-degree > d — do + 1. (11) 



i<i<d iej^ dQ iei 



(iv) Let zjc = ^2 ieI c Zi. Then 



e 2l+ "' +Zd = A(/)e 2/c mod terms of z-degree > d - d + 1. (12) 

Proof. Let Iq € ^>d - Then there exists ^ G Lj- Q such that £ € Lj- if and 
only if Lf o C Lj-, i.e., if and only if I C I . Thus (i) 44> (ii). 

Let us prove that (ii) (hi). We write Fi(zi) = 1 + Zigi(zi). We have 

l<i<d KC{l,...,d}ieK 
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Consider a monomial z^ ■ ■ ■ Zj d of total degree k\ + • • • + < d — d$. Let 
us denote its coefficient in the product ILeA" z i9i( z i) by ax- Let Kq be 
the set of indices such that fcj ^ 0. Then \Kq\ < d — d$. Moreover in 
the right hand side of (113D . our monomial appears only in the terms where 
K C i^ - Therefore the coefficient of z^ 1 ■ ■ ■ z d d in Oi<i<d Fi( z i) is equal to 
S.ft:c.ft: Furthermore, the coefficient of ^ 1 • • • z^ d in the right-hand-side 
of (fTtj) is equal to 

£ A(J) £ ajc = £ a* £ A(I). 

/ e j-d KCK ni c kck IeJ d 

— "0 —"0 

By condition (ii) we have 

^2 X(I) = 1 for every K C K - 

Thus we have proved that (ii) =>■ (iii). Next, (iv) is a particular case of (iii), 
so it remains only to prove that (iv) implies (ii). 

By expanding the exponentials in condition (iv), we obtain 

{z l + --- + z d f-^= £ A(/)(5» d ~ d °. 

Condition (ii) follows easily from this relation. □ 

4.2. Formula for intermediate sums. In preparation for the approxima- 
tion theorem, we need some notations and an expression for intermediate 
sums S L (s + c, A)(£). 

We have V = Lj © L/c For x GV we denote the components by 

X = Xj + Xjc 

Thus we identify the quotient V/Lj with Ljc and we denote the projected 
lattice by Ajc c ■ Note that Ljc n A C A/c , but the inclusion is strict in 
general. 

Example 22. Let v x = (1,0), v 2 = (1,2), 1= {2}. The projected lattice A/c 
on L/c = M Wl is Z^. See Figure H 

We denote by c/ the cone generated by the vectors Vj, for j & I and 
by b/ the parallelepiped b/ = ^2 i£ j[0,l[vi. Similarly we denote by C/c the 
cone generated by the vectors Vj, for j € I c and b/c = ^2 iGl c[0,l[vi. The 
projection of the cone c on V/Lj = Ljc identifies with c/c. Note that the 
generators Vi, i € I c , may be non-primitive for the projected lattice A/c, 
even if it is primitive for A, as we see in the previous example. We write 
s = si + sic. 

We first show that the intermediate generating function S Ll (s + c, A) 
decomposes as a product. 
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'iii i i i i i i i 

'ill I I I I I I I 

1 I I I I I I I I I I I 

• ••••• 

I I I I I I I I I I I 
I I I I I I I I I I I 



Figure 2. The projected lattice 

The function S(sic + c/c, A/c)(^) is a meromorphic function on the space 
(Ljc)*. The integral I(si + cj,Lj n A)(£) is a meromorphic function on the 
space (Li)*. We consider both as functions on V* through the decomposition 
V = Li® Lie 

Proposition 23. The intermediate sum for the full cone s + c breaks up 
into the product 

S L '(s + c, A)(0 = S{sjc + c/c, A /c )(0 I(si + CJ) Lj n A)(£). (14) 

Proof. The projection of the cone s + c into L/c is the cone s/c + c/c. For each 
xic G (s/c + c/c) n A/c, the slice (s + c) n (x/c + Lj) is the cone x/c + sj + c/. 
Let us compute the integral on the slice. 

/ e«*> dm L/nA (y). (15) 

J(s+c)n(x /c +Lj) 

We write y = xjc + s/ + X^ g / J/j^j- Then 

dm L/n A(y) = volL jn A(MlI d ^'' 
Hence (JT3J) is equal to 

e <^/c) e «,.r> vol L/ nA(b7)(-l) |/1 /T^T' 

We observe that only the first factor, e^' Xlc ' , depends on x/e. The sum of 
these factors over all x/c e (s/c + c/e) n A/e gives £(s/c + c/c, A/e)(£), and 
using formula ([I]) for the integral, we obtain (|14h . □ 

4.3. Approximation theorem. We can now state and prove the approx- 
imation theorem. 

Theorem 24 (Approximation by a patched generating function). Let c CV 

be a rational simplicial cone with edge generators Vi, . . . , Vd- Let s € Vq. Let 
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1 1->- A(I) 6e a patching function on J7>^ ■ Lor / € J7> d iei L/ 6e i/te linear 
span of{vi}i£i. Then we have 

S(s + c,A)(O^A\s + c,A)(0-= E A(7)^( S + c,A)(0 

mod terms of ^-degree > — c?o + 1- (16) 

We call the function A x (s + c, A)(£) on the right-hand side of (|16p the 
patched generating function of s + c (with respect to A). 

Proof of Theorem\24\ We write the vertex as s = Y^i s i v i- Let a = ^ OjUj G 
V. We apply (fTTI) to the functions 



F .( 2 .) = e (*iSi>i 



-Z; 



1 — e Zi 

and we substitute Zi = (£, i>j). We obtain 

11 1 _ e <£,vi) v ; 11 i_e«,»i) 

j=i ie^,, *g/ c 

— "0 

mod terms of ^-degree > d — do + 1. 

We multiply both sides first by e^' s ^; because this is analytic in £ and thus 
of non-negative ^-degree, the identity modulo terms of high ^-degree still 
holds true. Then we multiply by 1/ nf=i( — (£> u t))> which is homogeneous 
of degree — d in £. We obtain 



mod terms of ^-degree > —do + 1. (17) 

Now, we sum up equalities (fTT|) when a runs over the set (s + b) n A of 
integral points in the fundamental parallelepiped s + b of the affine cone 
s + c. On the left-hand side we obtain 

d 

V TT - 

ae(s+b)nA i=l 

By Lemma @ this is precisely S(s + c, A)(£). On the right-hand side, for 
each /, we have a sum over o € (s + b) n A of the function 



11 l_ e <€,«i) e 1J -(£,w 



which depends only on the projection ajc of a in the decomposition a = 
07- + ajc G Lj © L/c When a runs over (s + b) PI A, its projection ajc runs 
over (s/c + b/c) n A/c. Let us show that the fibers have the same number 
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of points, equal to vol/^nA^/). For a given a/c e (s/c + bjc) n A/c, let us 
compute the fiber 

{y € (s + b) n A : y/ c = a/c}. 

Fix a point a/ + a/c in this fiber. Then y = ajc + yj lies in the fiber if and 
only if yj — ai £ (s/ — ai + bj) D A. By Lemma Efii) , the cardinality of the 
fiber is equal to voIljciaO^i)- Thus we obtain 

S(s + c)(0 = 

E A ( J )^ + M(o n ri^ e<c ' s/> voi^nA(bz) n 

ltzJ >d 

mod terms of ^-degree > — do + 1. (18) 

By Proposition [23 and Lemmas [7] and EJ the term corresponding to an 
I G ^>d m the right-hand side of (PT8l) is precisely 

S Ll (s + c, A)(0 = 5(s Z c + c/c, A /c )(0 /( SJ + c/ , Lj n A)(£), 

which completes the proof. □ 

Remark 25. For do = 0, we obtain the poset J>§ of all subsets of {1, . . . , d}. 
The unique patching function on J7> is given by A(0) = 1 and X(I) = 
for all 1^0. T/ien £/ie approximation is trivial, i.e., S{s + c, A)(£) = 
A\s + c,A)(0. 

Example 26. Let c be the standard cone in K 2 , and do = 1- Thus consists 
of three subsets, {1},{2} and {1,2}. A patching function is given by A({i}) = 1 
and A({1,2}) = — 1. We consider the affine cone s + c with s — (— 5,— |). Let 
€= (6,6)- We have 

_ c -«>/2 e -Cl/2-? 2 /2 

/(«i + c {<} )(0 = — 7 , l(s + c)(0 = T7 . 

Si S1S2 

^ + - T^' S{s + cm = (i-^Ki-^) - 

The approximation theorem claims that 

1 1 -e"? 1 / 2 1 -e"« 2/2 e -« l/2 ^ 2/2 



(l_efi)(l-ee>) l-e« 2 a 1 - e& 6 66 

mod terms of ^-degree > 0. 



Indeed, the difference between the two sides is equal to 

~ 6/2 \/ 1 e" 



/ 1 e-^w 1 e"« 2 / 2 



which is analytic near 0. 
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4.4. An explicit patching function. Next we compute an explicit patch- 
ing function on J7>^ • It is related to the Mobius function of the poset J7> rfo , 
so we call it the Mobius patching function and denote it by AmoMus- We 
will denote the corresponding patched generating function A XM ° hius ( s _|_ C; a.) 
by ^> do (s + c,A). 

Proposition 27. For I G jt , , let 



>d 

I\-l 
d -l 



AMobius(/) = (-l) |/hd0 ^ 



Then XmoUus *s a patching function on ■ 

Proof. We prove that AmsMus satisfies Condition (iv) of Lemma [2TJ The 
trick is to write e z = 1 + t(e z — l)|t=i- Thus 

d d 

e zi+-+z d =ne 2 =n( i+ t{ ? zt - *)) 



1 i=l 



t=l 



Let us consider P(t) := nf=i( 1 + *(e^ i -1)) = Eg=o C q( z ) tq as a polynomial 
in the indeterminate t. As e Zi — 1 is a sum of terms of Zj-degree > 0, we 
have 

ko 

e zi+-+z d = y^^Cg(z) mod terms of z-degree > k + 1. (19) 
Next, we write 

d d 

pit) = n( i + ^ - !)) = ik* - o + tezi )- 
i i 

By expanding the product, we obtain 



\K\<q 

Summing up these coefficients for < q < k$ = d — do , we obtain 

ko I ko 



Ew = E E^r 1 * 1 



d~ \K\ 
- \K\ 



9=0 |K"|<feo \<3=l^l 

By substituting K = I c and d — q = m, we obtain 

ko 

Y,c q {z)= e /(i J i)^ c , 

9=0 \I\>d 

with 

m=do m=do 
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The truncated binomial sum f(j) is easy to compute, using the recursion 
relation ( j ) = + We obtain 

Thus, Am6Mus(-^) = f(\I\) satisfies Condition (ii) for a patching function. □ 



Remark 28. Proposition 27 can also be deduced from results in |16| . 



5. Computation of the patched generating function 

In this section, we show that if ko = d — do is fixed, the patched gen- 
erating function A>d Q (s + c, A) can be efficiently computed for a simplicial 
cone s + c. This will be a consequence of Barvinok's polynomial-time de- 
composition of cones in fixed dimension OUT]. We exhibit the dependence 
of the patched generating function on the vertex s explicitly as a "step func- 
tion" in two useful ways, using the "ceiling" function [•] and the "fractional 
part" function {•}, respectively. 

We will write the patched generating function using the following analytic 
function. 

Definition 29. Let 

oo n 

T(r,x) = e™-^- = - J> n (T)^, (20) 

n=0 

where B n (T) are the Bernoulli polynomials. 

We start with the following result. 

Theorem 30 (Short formula for S Ll {s + c,Z d )(£) for varying s). Fix a 
non-negative integer k$. There exists a polynomial time algorithm for the 
following problem. Given the following input: 

(11) a number d in unary encoding, 

(12) a simplicial cone c = c(«i, . . . ,Vd) C M. d , represented by the primitive 
vectors v 1 , . . . , S Z d in binary encoding, 

(13) a subspace Lj = lin(t>j : i £ /) C M, d of codimension ko, represented 
by an index set I C -[1, . . . , d} of cardinality do = d — ko, 

compute the following output in binary encoding: 

(01) a finite set T, 

(02) for every 7 in T, integers , rational vectors rjp and for 
i = 1 , . . . , d, where rj^ G 7L d for i G I c 
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such that for every s G Q d , we have the following equality of meromorphic 
functions of £: 

s L '{ S +t,z d m 

=E« (7) ^nf(^ ) ^}l,(^! 7, )) 

7er iei c 



nti<e, 



w 



= e«.'> £ a « II r( { - fa?" , «> } , <e, ^ » • . / (7) • (21b) 
7er ie/ c lli=i\4)^ / 

Of course, for / = we have L = {0}, and so we recover formulas for 
S(s + c,Z d )(£). If we set I = {l,...,d}, then L = R d , and we get formulas 
for /(s + c,Z d )(£). 

Remark 31. Consider the term corresponding to 7 G T in (|21a|) or (|21bD . 
^4s it toZZ follow from the proof the vector for i £ I is just the original 
vector Vi, and the collection i = l,...,d, forms a basis o/M d . Fur- 

thermore the vectors wp\ with i G I c , are in Ljc and form a basis of the 
projected lattice. The vectors rj^, i = 1, . . . ,d, are the dual (biorthogonal) 
vectors to the elements , j = 1, • • • ,d, i.e., (rji > w j ) = Thus we 
only need to compute the integers and the elements where i G I c . 

Remark 32. Consider the term corresponding to 7 G V in (|21b|) . As the 

vectors i G I c , form a basis of the projected lattice, we may identify 

V/(Lj+A) to ® ieI c[0,l[wp ] . Define 

* (7) = £{-<^>K (7) - 

iei c 

As the np^ for i G I c are integer vectors, and (ry| 7 ' \vj) = if j G I, 
we can think of s 1— )■ s^) as a linear map on the torus V/(A + Lj) with 

integer coefficients. The point s + s^ is in ® i6 j c Zi<;j © ® igJ Mvj, and 
formula (I21bjl reads also 



S^( S + a d )(0 = V« W e^ +s(7) > n - . (22) 



s (7) 

5tr n fc 7c(l-e<^"))n ie7 (e,^) 
Now we prove the theorem. 



Proof of Theorem \3(A Let us describe the algorithm along the proof. Let 
A = Z . By Proposition! 



S Ll {s + c, A)(0 = S( Sl c + c/c, A /o )(0 /(sj + c Jt Lj n A)(0- (23) 
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We first discuss I(sj + Cj, Lj fl A). We have 

I(si + cj, Lj n A)(0 = vol L/ nA(b/) TT -jf^-. (24) 

Using linear functionals rji € Q d , i E I (the coordinate functions with re- 
spect to the basis write sj = ^2 i€ j{r/i, s)uf. The r/j can be read off in 
polynomial time from the inverse of the matrix whose columns are t>i, . . . , 
Then e^' Sl ^ takes the form 

e^)=nexp((r ?i , S )(e,^)). (25) 

Now, to handle the factor S(sic + c/c,Ajc), note that C/c c Ljc is a fco- 
dimensional cone. By using a Hermite normal form computation, which is 
polynomial time [29], we can compute a linear change of variables which 
replaces the projected lattice Ajc on Ljc by Then, using Barvinok's 

decomposition [9], we decompose it into a family of cones which are uni- 
modular, 

[c/c] = Cmfc/™''] (modulo cones containing lines), (26) 

mGM 

where e m € {±1}- As /co is fixed, this decomposition can be done by a 
polynomial time algorithm. This step is of course crucial with respect to 
the efficiency of the whole algorithm. 

Changing again notations, we now denote by c = c({u;j}i e jc) one of these 
unimodular cones t™^ C -L/c, with primitive generators Wi, and also write 
e = e m . We remark that the vectors Wi, i € I c , generate the projected lattice 
on Ljc. Using linear functionals r\i 6 Q d , i € I c , write sjc = Y^ieleiWii 3 ) w i' 
Actually, we have rji £ Z d . By letting Wi = V{ for the other indices i S /, we 
can write 

d 

s = } ](r)j,s)vi + } j {rn,s)w i = } j {r)i,s)w i . (27) 

Let Sj c be the unique lattice point in the fundamental parallelepiped of 
the cone sjc + c. We have 

^c = X)r<»K.*>M- (28) 
Using this, we obtain the generating function from Lemma [9] as 

Thus finally, using (I24D we have the meromorphic function 
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Then (i29j) is now written 

« n Tfl-fa,*)!, • n ex p(^' s )(^^)) • 77^-7; — r. (3°) 



where a collects the multiplicative constants in (|29|) . Collecting these terms 
gives the desired short formula (|21ap . 

To derive the second form, we note \(rji,s)~\ = (rji,s) + {— (r)i,s)}, so we 
can write 

T(\(th, s)),(C, Wi)) = T({-{ m , s}}, <£, Wi))exp({rH, s)<£, Wi )). (31) 
Thus the term (|30|) can be written as 

« II r({-<i», *>}, (£, vh)) • e«' s > • , \ - , (32) 

using (|27p . Collecting these terms gives the short formula (|21b[) . □ 

Example 33. Let us give an example of the output of our algorithm, in a small 
example. Consider the 3-dimensional cone with rays given by the vectors (1,1,1), 
(1,-1,0), (1,1,0). This cone is not unimodular. We consider the affine cone 
s+C. Our algorithm described in Theorem l30l computes any intermediate generating 
function S L (s + c, 1?) when L is a linear span of a face of c. For L = {0} (indexed by 
the empty set /), we obtain the meromorphic function S(s + c, Z 3 )(£) (the discrete 
generating function of the cone s + c). Here S(s + c, Z 3 )(£) depends of s = (s%, S2, S3) 
and £ = (fi: & i &) and is given by: 

exp (si£i + s 2 £ 2 + S3£3) 

T ({-s 3 + s 2 } , -£1 - £2) T {{-s x + s 2 } , £ x ) T ({-s 3 } , £1 + £ 2 + £3) 



+ 



(-£i-£ 2 )£i (£i+£ 2 +£ 3 ) 

T ({-.s 3 + s 2 | , a - £2) T ({2 .S3 - s 2 - si} , £1) T ({-s 3 | , £1 + £2 + £3) 
(£i-£ 2 )£i (£i+£ 2 + £ 3 ) 



If L = Mvi is the subspace of dimension 1 generated by the edge V\ = (1,1,1) of 
the cone c (so that L is indexed by the subset / = {1} of {1,2, 3}), the intermediate 
generating function S L (s + c,Z 3 ) is given by: 



exp(si£i + s 2 £ 2 + s 3 £ 3 ) 



T ({-s 3 + s 2 } , -£1 - £ 2 ) T ({-si + s 2 | ,£1) 

(-£i-£ 2 )£i(-£i-£2-£ 3 ) 
T ({-s 3 + s 2 | , £1 - £2) T ({2 s 3 - s 2 - si} , £Q 
(£i-£ 2 )£i(-£i-£2-£ 3 ) 



Remark 34. When k$ = d — is fixed, the set J^ d has a polynomially 
bounded cardinality, and it can be enumerated by a straightforward algorithm 
along with the evaluation of the patching function AmsWus • Thus we can also 
compute A>d { s + c > A)(£) in the same form (|21aj) or (|21b|) in polynomial 
time. 
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6. Computation of Ehrhart quasi-polynomials 

We now apply the approximation of the generating functions of the cones 
at vertices to the computation of the highest coefficients for a weighted 
Ehrhart quasi-polynomial. We first discuss the case when the weight is a 
power of a linear form. 

Theorem 35. Let p be a simple polytope and let V(p) denote the set of its 

vertices. For each vertex s E V(p), let c s be the tangent cone of s, and let 
q s G N be a positive integer such that q s s € A. Fix a linear form £ £ V* 
and M a nonnegative integer. Fix < k$ < d and let do = max{d — ko,0}. 
Then the Ehrhart quasi-polynomial 

E(p,£,M;n) = £ ^ 
xgnpnA 

coincides in degree > M + d — ko with the following quasi-polynomial 
k ° It „\M+d-k 

£ E (H,J M+d " fc ,%! d _ kV ({»}«.* + c„ A) M+fc] (C), (33a) 

fc=0s€V(p) V '' 

evaluated at £ = I, which can also be written as 

k=o se v{p) y ' 11 

(33b) 

evaluated at £ = £. 

In the following, we will use the second form (|33b|) . 

Remark 36. The sum (|33h depends polynomially on i. However, for an 
individual vertex s, the functions 
It g \M+d-k 

* * (M + d-fc)! A> ~ d ° {{nUS + Cs ' A) 

and 

It „\M+d-k . N 
5- (e-'«"'->^ o({ „ } „ s + c.,A)(o) M+t| 

are meromorphic functions, which are not defined if £ is singular. Thus in 
the algorithm we use a deformation procedure. 

Proof of Theorem [35l The sum Y^xenpnA * s * ne term of ^-degree M 

in 

S(np)(0= E S(ns + c s )(0- 
sev(p) 

Fix a vertex s. We write n = [n\ qa + {n} qs . As L n J<? s s is a lattice point, we 
have 

S(ns + c s )(0 =eW«.«.»>5({n}„ S + c a )(0. (34) 
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Consider S(ns + C s )nvfi(£) as a quasi-polynomial in n. By it coincides 
in degree > M + d — /co with 

k o it \M+d-k 

BN,j AW ( ^^_ fc) , s(w 9sS + c s ) M+fe] (o. 

Now, for < k < /cq, we have 

£(W?.« + ^)[-d+fe](0 = A >d {{n} qs s + C a )[-d+.k](0- 
By specializing on ^ = £, we obtain the claim in the form of equation (|33ap . 
To obtain the second claim in the form of (|33b|) . we write 

S(ns + c s )(0=e n ^ s \e-^ n ^S({n} qs s + c s )(0). (35) 

Again, by expanding we obtain that the quasi-polynomial S(ns + c s )[m](C) 
coincides in degree > M + d — k$ with 



Since e~^' s ^ n ^ qa is analytic in £, we have for < k < ko that 

L-(M{n }qsS{{nUs + Csm \ 

\ / [— d+k\ 

= L-(ZMn} qsA>do{{n}g3S + Csm ) (36) 
\ / [— a+rej 

Again, by specializing on £ = £, we obtain the claim in the form of 
equation (|33bj) . □ 

We now derive the coefficients of the weighted Ehrhart polynomial as 
short closed formulas that are "step polynomials" (cf. [34J). These can then 
be evaluated efficiently, providing a corollary (Theorem l42|) in the same form 
as Barvinok's theorem in 1101. 



Theorem 37. For every fixed number ko E N, there exists a polynomial-time 

algorithm for the following problem. 

Input: 

(11) a number d E N in unary encoding, with d > ko, 

(12) a finite index set V, 

(13) a simple poly tope p, given by its vertices, rational vectors Sj E Q d for 
j E V in binary encoding, 

(14) a rational vector I E Q d in binary encoding, 

(15) a number M E N in unary encoding. 

Output, in binary encoding, 

(01) an index set T, 

(02) polynomials p ,m E • • • , r k ] an d integer numbers Q' E Z, 
qj' m E N for 7 E T and m = M+d — ko,... ,M+d and i = 1, . . . , ko, 
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such that the Ehrhart quasi-polynomial 

ip v m M+d 

E(p,£,M;n) = £ = ^ E m (p,£,M;{n} q )n m 

xGnpnA m=0 

agrees in n-degree > M + d — ko with the quasi-polynomial 

M+d 

E E / 7,m ({crn}^, . . . , {«™»><0 » m - 

76r m=M+d—ko 

Remark 38. For d < fco, i/ie algorithm actually computes the complete 
Ehrhart quasi-polynomial, i.e., the coefficient functions E m (p,£,M;{n} q ) 
for m = 0, . . . , M+d. The key point of our method, however, is to handle the 
case where d > ko; then the non-trivial efficiently computable approximations 
come into play. 

Remark 39. The specific form of the quasi-polynomial given by the theorem 
gives a more precise period % for the individual terms, rather than a period 
q s that is determined by the vertex. The qi will always be divisors of q s . Due 
to the projections into lattices in small dimension < ko, these periods can be 
much smaller than q s . In particular, the highest- degree coefficient EM+d of 
course is a constant. 

We will use the following lemma. 

Lemma 40 (Lemma 4 of [2]). For every fixed number D G N, there exists 
a polynomial time algorithm for the following problem. 

Input: a number M in unary encoding, a sequence of k polynomials Pj G 
Q[Xi, . . . , Xd] of total degree at most M, in dense monomial representation. 
Output: the product P\- ■ ■ truncated at degree M. 

We can now prove the theorem. 

Proof of Theorem \37\ Because the polytope p is simple, we can use the 
primal-dual algorithm by Bremner, Fukuda, and Marzetta |17} Corollary 1], 
to compute the inequality description (H-description) from the given V- 
description in polynomial time. From the double description, we can com- 
pute in polynomial time the description of the tangent cones c Sj . for j G V by 
the primitive vectors v Sj7 i, ... , v Sji d & %> d such that c Sj = c(v Sjt i, . . . , v Sj ^)- 
We now use formula (|33b[) of Theorem [35l which gives (with do = d — ko) 

E m (p,^M; {n} q ) 

- E ^r( e - <w *-* > ^.(("i^ + c - A )«)) M+fc| < 37 > 

sev(p) 

for m = M+d—k, when m > M+d— ko. We compute this separately for each 
k = 0, . . . , ko, that is, m = M + d — ko, ■ ■ ■ , M + d. Let s + c s be one of these 
cones. By the algorithm of Theoreml30l and Remark [341 we compute the data 
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describing the parametric short formula (|21bj) for A>d {{ n }q a s + c sjA)(£). 
We then consider one of the summands of 

M^( e -«,M,..>^ ({ „ }jiS + c „ A)K) ) [ _ j+t] 

at a time. Here e~^^ n ^is s ) an d the term e^'^ n ^ qaS ^ from (]21bp cancel, and 
thus each summand takes the form 

( I |^) ( nnn(n)x^o))j fE ^) (-) 

where 

Ti(n) :={-<»*, s){n} ff .} for i € I c . (39) 

Let E N be the smallest positive integer such that qi{—r]i, s) E Z. Then <ft 
is a divisor of the number q s associated with the vertex s, because r\i E Z d . 
Then 

T *( n ) = il&MgJ^ with & = Qi(-Vus) EZ. 
Since % is a divisor of g s , this simplifies to 

*(») = itt**}., ( 4 o) 

where of course Q can be reduced modulo % as well because n is assumed 
to be an integer. We now treat n := {0n} qi E N as symbolic variables. 

In order to evaluate (j38[) at £ = £, we note the following. The first factor 
is holomorphic in £ and homogeneous of ^-degree m = M + d — k, the 
second factor is holomorphic in £ and homogeneous of degree k, and the 
third factor is homogeneous of ^-degree —d. If (£,Wi) = for some i, we 
cannot just substitute £ = £ in the formula. Instead we use a perturbation. 
In polynomial time, we can compute a rational vector £' E Q d such that 
(£',Wi) 7^ for all vectors u;, with (£,Wi) = 0. It is important that we 
choose the same vector once and for all computations with all cones and 
summands. 

We then set £ = t(£ + e£'), where t and e are treated as symbolic variables. 
Here the exponent of the variable t keeps track of the ^-grading. We then do 
computations with truncated series in Q[r^ : i E I c ] [t , e^ 1 ] . We note that 
this is a polynomial ring in a constant number of variables only, because \I C \ 
is bounded above by the constant &q. Thus Lemmal40l gives us a polynomial- 
time algorithm for multiplying the series. Then (|38p can be written as: 

(£ + e£',s} m (t~I rr(-..r~-\ Urn i -e'\ „.. \^ 1 t M ~ k 



, TT T(ri(n),(t(£ + e£'),Wi))) ■ —, 

(41) 

where the subscript [k] now means to take the term of t-degree k. In the 
end we are interested in the coefficient of the term t M e°. 

Expanding the factors of (|4ip gives the following contributions, all of 
which can be written down in polynomial time. First of all, the rational 
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terms (£ + e£',Wi) 1 give the following contribution. If (£,Wi) = 0, we 
simply get 

e~\ (42) 



(£ + e£', Wi ) {£', Wi ) 
If {£, Wi) / 0, we get the geometric series in e 



(£ + e£>,Wi) (£,Wi)^ \ (£,w. 

The first and second terms in (j41|) are holomorphic, thus the only negative 
degrees in e come from the rational terms ()42p . Let U be the number of 
vectors W{ that are orthogonal to £; then e~ u is the lowest negative degree. 
Note that U < d. Since we wish to find the term of e-degree 0, we can 
truncate all series after e-degree U : 

i i (£',wi 



(£ + e£', Wi ) (£,Wi)^ \ (£,w. 



We expand the first factor of (|4ip as follows. 



m! 

M=0 



S ( ) S ) m ""^' + °^ eC/ )- ( 44 ) 



Now we consider the holomorphic terms 
T( n (n),{t(£ + e£'), Wi )) 



oo 1 




7 

, --0 J 

/ min{j,C/} . 

^ ' [Lw,)'-" [('.«<,)"<" | I' 

(45) 

The Bernoulli polynomials Bj(ji) of degree j < ko that appear in this for- 
mula can be efficiently expanded in polynomial time using recursion formu- 
las. We remark that the variables appear with a degree that is at most 
that of t. Using Lemma HUl we multiply the truncated series (|45p for i € I c 
in Q[rj : i € -f c ][£][e], truncating in each step after t k ° and e u . We thus 
obtain the second factor of (|4"T|) . 

T^nWt^ + eO,^})) for all fc = 0,...,fco, (46) 

truncated series in Q[rj : i E i~ c ][e]- 
Then we multiply the truncated series ((121) . (jl3|) . (116]) . and (jH|) in poly- 
nomial time, truncating in each step after e u , using Lemma HU1 In the end, 
we read out the coefficient of e° as a polynomial in Q[rj : i 6 / c ]. Then we 



(n 
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substitute for r^. Collecting these terms gives the formula for the Ehrhart 



Example 41. Let us give a small example of the output of our algorithm for 
E m (p, I, M, {n} q ), when p is the simplex in R 5 with vertices: 



(0,0,0,0,0), (i,0,0,0,0), (0,i, 0,0,0), (0,0, £,0,0), (0,0,0,1,0), (0,0, 0,0, J). 



We consider the linear form I on M 5 given by the scalar product with (1,1,1,1,1). 

If M = 0, the coefficients of E m (p,£, M — 0; {n} q ) are just the coefficients of the 
unweighted Ehrhart quasi-polynomial S(np, 1). We obtain 



Now if M = 1, all integral points (x\, x 2 , £3, £4, x 5 ) are weighted with the function 
/i(x) = x\ + x 2 + x 3 + X4 + £5, and we obtain 



We can remark that although g = 6 is the smallest integer such that qp is a 
lattice polytope, only periodic functions of n mod 2 enter in the top three Ehrhart 
coefficients. This is indeed conform to the known periodicity properties of the 
Ehrhart coefficients. 

As a corollary, simply by evaluating the step polynomials, we obtain the 
following result, which directly extends the complexity result from Barvi- 
nok's paper to the weighted case. 

Theorem 42 (Evaluation of the Ehrhart coefficients for a given dilation 
class {n} q ). For every fixed number ko £ N, there exists a polynomial-time 
algorithm for the following problem. 
Input: 

(11) a number d € N in unary encoding, with d > ko, 

(12) a finite index set V, 

(13) a simple polytope p, given by its vertices, rational vectors Sj € Q d for 
j € V in binary encoding, 

(14) a rational vector £ G Q d in binary encoding, 

(15) a number M € N in unary encoding, 
(Iq) a number n in binary encoding, 

Output, in binary encoding, 

(01) a positive integer q € N such that qp is a lattice polytope and 

(0 2 ) the numbers E m (p, £, M; {n} q ) for m = M + d - k , . . . , M + d. 

Remark 43. A direct algorithm for computing E m (p, £, M; {n} q ) for just 
one dilation class {n} q could of course use the values r,- L = {Cin} qi € Z rather 



coefficient E m (p, £, M; {n} q ). 



□ 
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than symbolic variables T{ and would therefore only need to do calculations 
with truncated series in the two-variable ring Q^^e ]. 

Via the decomposition of polynomials into powers of linear forms, which 
is, as discussed in [2j, polynomial-time under suitable hypotheses, we obtain 
the following corollary. 

Corollary 44. For every fixed number ko € N, there exist polynomial-time 

algorithms for the following problems. 

Input: 

(11) a number d € N in unary encoding, with d > ko, 

(12) a simple rational polytope p C M. d , given by its vertices in binary 
encoding, 

(13) a number M in unary encoding, 

(14) a polynomial h of degree < M which is given either as 

(a) a power of a linear form, or 

(b) a sparse polynomial where each monomial only depends on a fixed 
number of variables, or 

(c) a sparse polynomial of fixed total degree, 

(15) a number n in binary encoding, 

Output, in binary encoding, 

(01) a positive integer q £ N such that qp is a lattice polytope and 

(02) the numbers E m (p, h; {n} q ) for m = M + d — ko, . . . , M + d. 

7. Experiments 

We implemented the algorithms in Maple, for the unweighted case and 
assuming that the input were lattice simplices of full dimension (in this case 
the quasi-polynomial becomes a polynomial). This assumption was made 
for simplicity of output in the calculation and because available software to 
verify the results (e.g., LattE macchiato [31]) cannot compute with weights. 
In addition, already the problem of computing Ehrhart polynomials for lat- 
tice simplices has received attention by many researchers and it is non-trivial 
(see e.g., the references in |21j). After checking simple low-dimensional ex- 
amples by hand, we set up automatic scripts for generating random tests. 
The simplices generated had vertex coordinates drawn uniformly at random 
from {—99, . . . ,99}. We timed the speed of the procedure to compute the 
top three Ehrhart coefficients in 50 random simplices per dimension and 
recorded the average time of computation. We compared with the compu- 
tation of the full Ehrhart polynomials using the state-of-the-art algorithms 
implemented in LattE macchiato [31] ; see Table [TJ 

In the table, Dual refers to an implementation of Barvinok's decomposi- 
tion of the duals of the tangent cones into unimodular cones, as implemented 
first in LattE |25j , and which is still the default method in LattE macchiato^ 
Primal refers to a primal variant of Barvinok's decomposition described 



The LattE macchiato command is count — ehrhart-polynomial. 
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Table 1. Computation times for Ehrhart polynomials of 
random lattice simplices 



Average runtime (CPU seconds) 





Full (LattE macchiato) 


Top 3 
(new code) 


Dimension 


Dual 


Primal Primaliooo 


3 


0.16 


0.10 0.04 


1.12 


4 


28.00 


4.68 0.28 


4.31 


5 




317.5 5.8 


13.4 


6 




198.0 


37.4 


7 






103 


8 






294 


9 






393 


10 






1179 


11 






1681 



in |30j : it is more efficient for these examples because the determinants of 
the dual cones are much larger^ We remark that our implementation of 
the new algorithm in Maple also uses a primal variant of Barvinok's de- 
composition to unimodular cones, which was introduced in [20] . Thus the 
new code should be compared to the runtimes listed in column Primal. Fi- 
nally, Primaliooo refers to a variant in which Barvinok's decomposition is 
stopped when a cone has a determinant at most 1000; then the points in 
the fundamental parallelepipeds are enumeratedJl 

All computations were stopped if unfinished after 30 minutes, thus the 
table ends at dimension 11 because all randomly generated examples we tried 
in dimension 12 took more than 30 minutes of calculation. The computation 
times are given in CPU seconds on a computer with AMD Opteron 880 
processors running at 2.4 GHz. 

In conclusion, the experiments indicate that the algorithms presented here 
can lead to dramatic improvements upon the computation of full Ehrhart 
polynomials. The fact that, for very low dimensions, the implementation 
is slower than LattE macchiato, is explained by the choice of Maple as an 
implementation language. Maple is an interpreted system, which is much 
slower than C++, the implementation language of LattE macchiato. We 
expect that the speedups of Primaliooo compared to Primal, which were 



The command is count — ehrhart-polynomial — irrational -primal. 

n 

The command used is count — ehrhart-polynomial — irrational -primal 
— maxdet=1000. 
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first documented in [30], will also be obtained in a refined implementation 
of our new algorithms. 

The implementation is available at [3]. 
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